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Abstract Most dusty winds are described by a set of similarity functions of a single inde- 
pendent variable that can be chosen as rv, the overall optical depth at visual. The 
self-similarity implies general scaling relations among the system parameters, in 
agreement with observations. Dust drift through the gas has a major impact on 
the structure of most winds. 
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1. Introduction 

The complete description of a dusty wind should start with a full dynamic 
atmosphere model and incorporate the processes that initiate the outflow and set 
the value of M. These processes are yet to be identified with certainty, the most 
promising are stellar pulsation (e.g. Bowen 1989) and radiation pressure on the 
water molecules (e.g. Elitzur, Brown & Johnson 1989). Proper description of 
these processes should be followed by grain formation and growth, and sub- 
sequent wind dynamics. Two ambitious program attempting to incorporate as 
many aspects of this formidable task as possible have been conducted over the 
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past few years by groups at Berlin (see J.M. Martin, these proceedings) and Vi- 
enna (see Dorfi et al 2001). While much has been accomplished, the complexity 
of this undertaking necessitates simplifications such as a pulsating boundary. In 
spite of continuous progress, detailed understanding of atmospheric dynamics 
and grain formation is still far from complete. 

Fortunately, the full problem splits naturally to two parts, as recognized 
long ago by Goldreich & Scoville (1976). Once radiation pressure on the 
dust exceeds all other forces, the rapid acceleration to supersonic velocities 
decouples the outflow from the earlier phases — the supersonic phase would 
be exactly the same in two different outflows if they have the same mass-loss 
rate and grain properties even if the grains were produced by entirely different 
processes. Furthermore, these stages are controlled by processes that are much 
less dependent on detailed micro-physics, and are reasonably well understood. 
And since most observations probe only the supersonic phase, models devoted 
exclusively to this stage should reproduce the observable results while avoiding 
the pitfalls and uncertainties of dust formation and the wind initiation. 

2. Dusty Winds: Problem Setup and Solution 

A steady-state dusty wind is controlled by three forces. First, radiation 
pressure drives the outflow. When ry > 1, the radiation is continually reddened 
and its coupling to the dust is degraded; this weakening of the force can be 
described by a radial profile (/>(r). Radiation pressure acts mostly on the dust 
grains while the wind mass is mostly in the gas particles, and coUisional coupling 
between these two components introduces the second force, drag. Oilman 
(1972) has shown that the dust-gas relative velocity reaches steady state within 
a distance much shorter than all other scales. The drag force can then be 
eliminated, because the dust fully mediates to the gas the radiation pressure 
force. But the drift still has an important effect, since the radiative force per 
unit mass is proportional to nj:/n, and separate mass conservation for the dust 
and gas implies n^lno^vjv^} Finally, the gravitational pull of the central star 
opposes the expansion and must be overcome in any outflow. 

Accounting for the three forces, the equation of motion in terms of dimen- 
sionless velocity w and radial distance y is (EUtzur & Ivezic 2001) 



The constants P and F are the ratios of radiation pressure to the drag and gravity, 
respectively, C,=vlv^'\^ the drift profile and (\) is the reddening profile, obtained 



Even with prompt dust formation and no further grain growth or destruction, the dust abundance varies in 
the shell because of the difference between the dust and gas velocities. 




(1) 
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from the separate, but coupled, equation of radiative transfer. The outflow 
problem is fully described by this simple mathematical equation with only 
three free parameters — the constants P and F and the initial value uj{y = 1). 
There are no other free parameters. The wind structure will be exactly the same 
in systems that have different mass (M), luminosity (L) or mass loss rate (M) 
when these properties combine to produce the same P and F. 

The complete solution of equation 1 is presented in Elitzur & Ivezic (2001), 
including the analytic solution for optically thin winds (ry < 1). Here we 
present observational impUcations. 

2.1 Physical Domain 

The relation Mv^ < L/c has often been used as a physical bound on 
radiatively driven winds, even though the mistake in this application when 
TV > 1 has been pointed out repeatedly (e.g. Ivezic & Elitzur 1995). Instead, 
the proper form of momentum conservation is Mv^o = TpL / c, where rp is the 
flux-averaged optical depth, and it yields 

Mv^ = - {Q*/Qv) TV (1 + Tv)-°-^' . (2) 
c 

Here Qy is the dust efficiency coefficient at visual and is its Planck-average 
with the stellar temperature. Typical values of QJQ\ are 0. 1 for silicate and 
0.25 for carbon dust. Since the ratio Mvqo /{L/c) can exceed unity and increase 
without bound, momentum conservation does not constrain the wind. Instead, 
the constraints come from force considerations — the outward force must exceed 
gravity (positive derivative in eq. I), leading to the two phase space boundaries 
shown in figure I. The liftoff bound, shown in the left panel, sets a lower hmit on 
M, proportional to M^/L; below this minimal M the grains are ejected without 
dragging the gas because the density is too low for efficient dust-gas coupling. 
As a condition on the wind origin, this bound is the most uncertain part of the 
solution. The result M^iR oc /Lis reasonably secure (a similar relation was 
noted by Habing et al 1994), but the proportionality constant involves a fudge 
factor. This factor can only be determined from a more complete formulation 
that handles properly grain growth. The bound shown in the right panel reflects 
the weakening of radiative coupUng in optically thick winds because of the 
radiation reddening, setting an upper limit on M. 

2.2 Scaling 

Most stars are located well inside the allowed region of phase space, as is 
evident from figure 1. The dependence on F and w{l) disappears rapidly with 
distance from the boundaries, therefore most dusty winds are described by a 
set of similarity functions of just a single independent variable P, which is 
equivalent to ry. The self-similarity implies general scaling relations. Systems 



4 




Figure 1. Winds powered by radiation pressure on dust should fall above the liftoff bounds 
shown in the left panel and below the reddening bounds in the right one. L4 is the stellar 
luminosity in 10'* Lq (Ivezic & Elitzur 2002). 

with the same combination M jL^I^ necessarily have also the same Mvoq/L 
because both are uniquely related to ry- 

The velocity profile for all winds can be summarized with the simple expres- 
sion 

1_M' where |^ = i when tv < 1 
^ J , wnere | ^ ^ q 4 ^^^^ > ^ 

and where = ri[l — {vi/vooY^^], with vi the wind velocity at its starting 
radius ri. The profile shape of v/v^ is nearly the same for all winds, a weak 
dependence on ry enters only in the power k. The small-rv value A; = | reflects 
the effect of the drift, which dominates the dynamics in that regime, the switch 
to a more moderate profile at large ry is caused by reddening. 

The dimensionless solution involves only properties of individual grains. 
The dust abundance is irrelevant and affects only the velocity scale v^o. When 
only the radiation pressure force is taken into account, increases with L and 
is independent of M (Goldrecih & Scoville 1976). Drift effects change this 
behavior fundamentally, producing instead 

^;^ = ^M(l + Tv)-'■^ (4) 

where the proportionality coefficient A contains properties of the grain material 
and n^a^jn. Therefore, outflows with ry < 1 should obey oc M if they 
have the same dust properties. Remarkably, even though the wind is driven by 
radiation pressure, its velocity is independent of luminosity. The dependence 

on L enters only in optically thick winds because ry = A^/^{Qy/Q^,)M^^^/ L. 

The behavior implied by eq. 4 had been noticed by Young (1995) from 36 
nearby Mira variables with low mass-loss rates. Young's data is presented in 



The Structure of Winds in AGB Stars 



5 



36 stars from Young (1995) 
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Figure 2. Histograms of vi,/M for oxygen- and carbon-rich stars with rv < 1 outflows. 
In each data set, both Voo (in km s"'^) and M (in 10~® Mq yr~^) were determined from CO 
observations. A vertical line marks the mean of each distribution. The corresponding coefficient 
A (see eq. 4) is listed together with the value of (T22 it implies (Ivezic & Elitzur 2002). 



figure 2 as a histogram of v^/M togetlier witii a similar histogram for C-rich 
stars with small optical depths (data from Olofsson et al 1993). Each histogram 
shows a pronounced peak — as implied by eq. 4 when the dust properties do 
not vary among members of the set. The value of ndad/n = W~^^a22 cm~^ 
determined for each set from its average A and properties of silicate and carbon 
grains is listed in the figure. Remarkably, the two samples of different types of 
stars produce the same a22- In spite of the large differences in atmospheric and 
grain properties between O- and C-rich stars, the fraction of material channelled 
into dust is such that a22 turns out to be the same in both. 

Drift and reddening play a major role in the dynamics. As M decreases, 
the dust and gas decouple and the velocity decreases too. As M increases 
and the wind becomes optically thick, reddening degrades the efficiency of 
the radiation pressure force and the velocity again decreases. In between, the 
velocity reaches maximum fmax — 20L\^'^ km s^^, obtained at M{vma.x) — 

2x10^^X4'^^ Mq yr^^. The proportionality coefficients are slightly different 
for carbon and silicate dust. 

2.3 Individual Sources 

Since ry controls both the dynamics and the radiative transfer, in any given 
source the spectral energy distribution (SED) can be fitted with only one free 
parameter and the corresponding solution should provide a fit also for the radial 
velocity profile. Figure 3 shows the SED fit for W Hya. A byproduct of the 
solution that produced this SED is the shape of the dimensionless velocity 
profile. There is no freedom to modify this profile, and all the data points (from 
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Figure 3. Modeling the wind in W Hya (Zubko & Elitzur 2000) with outflow calculations by 
the code DUSTY (Ivezic et al 1999). Only one free parameter was varied, rv, determined to 
be 0.83 from the SED fit in the left panel. Right: Model prediction for the gas velocity profile. 
Data points show masers and CO thermal emission from the wind. The SiO maser emission 
originates in the extended atmosphere and is not part of the outflow. 

observations of CO thermal emission and various masers) indeed follow it when 
scaled to dimensionless values. Figure 3 shows the results with actual physical 
variables, setting the dimensional scales for r and v from the Hipparcos distance 
of 1 15 pc and the terminal velocity of 8 km s~^ from CO emission. 

Significantly, neither L, M or dust-to-gas ratio were input parameters of the 
model. Instead, these quantities are derived afterwards when the SED fitting 
results are supplemented by the distance and velocity scales. The bolometric 
flux is a direct product of the fitting process and together with the distance to 
the source determines the overall luminosity. And once L, v and tv are known, 
M and the dust abundance are calculated from general self-similarity relations. 

3. Water Thermal Emission 

Water is an abundant component of the outflow in oxygen-rich stars and the 
dominant coolant at radial distances up to ~ 10^^ cm. But because of atmo- 
spheric absorption, the discovery of its thermal emission had to wait for ISO. 
With the wind properties in W Hya determined from SED analysis, Zubko & 
Elitzur (2000) successfully fitted all 21 water Unes observed by ISO using H2O 
abundance as the only free parameter. Harwitt & Bergin (2002) subsequently 
found the same line ratios in VY CMa. This highly luminous supergiant has L 
= 5x10^ Lq and M ~ 2xl0~^ Mq yr~^ while in W Hya the corresponding 
values are only 1.1x10^ Lq and 2.3xl0~^ Mq yr~^. In spite of these large 
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Figure 4. Top Row: Contour plots of the observed CIT3 images at 1.24 /xm (J-band), 1.65 
/um (H-band), and 2.12 /im (K-band) from left to right (Hofmann et al 2001). Bottom: Model 
results of Vinkovic et al, 2002. 



differences, scaling the model results for W Hya by an overall factor of 10 
produces adequate fit for the H2O lines in VY Cma too. 

Thanks to the multitude of observed transitions, in H2O line analysis the 
abundance is derived rather than entered as input. This is a significant advantage 
over the common determination of M from CO observations, wliich involves 
only one transition and must assume an abundance beforehand. 

4. Deviations from Spherical Symmetry 

Many AGB shells display spherical symmetry, but there are numerous ex- 
ceptions. A peculiar case is the source CITS, one of the most extreme infrared 
AGB objects. Observations by Hofmann et al (2001) show a 1.24;um image 
elongated along a symmetry axis (fig. 4, top left panel). The elongation disap- 
pears at slightly longer wavelengths, and the image becomes almost perfectly 
circular at lAljim. This peculiar behavior has been successfully explained by 
Vinkovic et al (2002) with the aid of a newly developed 2D radiative transfer 
code, LELUYA (see www.leluya.org). In the model geometry, the wind spher- 
ical symmetry is broken by two bi-polar cones and the observing viewpoint is 
inclined to the axis. The density distributions in the conical regions and the 
bulk of the outflow are deduced directly from the observed brightness of the 
J-band image. The model results, shown in the lower panels of fig. 4, are hardly 
distinguishable from the observations. 
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The wavelength- variation of the images reflects interplay between scattering 
and emission. Scattered Ught dominates at A < 1 .5 ^m and the images trace 
the large-scale contours of the density distribution. This explains the elongated 
shape of the J-image. At longer wavelengths emission takes over and the images 
are affected by the shape of the dust isotherms, becoming circular due to the 
central heating of the dust. 

5. Conclusions 

The "standard model" seems to be working satisfactorily. Detailed fits are 
successfully obtained with a single free parameter, ry, as predicted. Dust drift 
is a major ingredient of the dynamics and the reason for the peculiar correlation 
discovered by Young (1995). 

Two fundamental issues remain open and can be addressed only by more 
complete models that incorporate the wind origin (such as those of the BerUn and 
Vienna groups). The first is a more reliable determinations of the phase space 
lower boundary shown in fig. 1. What is the minimal M? The second involves 
the dust abundance, which sets the scale of terminal velocities in AGB outflows 
to ~ 10 km s~^; were that abundance 100 times higher, typical velocities would 
be ~ 100 km s~^ instead. What is so special about nacrd/n ~ 10^^'^ cm^^ and 
why is it the same in both C- and O-rich stars? It is entirely possible that the 
answers to these two open issues are related. 
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